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ABSTRACT 


A numerical procedure was developed for minimizing the structural 
mass of an aircraft structure which must have a specified minimum 
flutter velocity or divergence velocity. During the optimization 
process the arrangement of the structural members remains fixed, while 
the stiffness parameters of the structure are varied. 


I. Recurrence Relation 


Let w be the total weight of an aircraft structure and let ^ be 
an aeroelastic response parameter such as the flutter velocity or 
the divergence velocity. It is desired to minimize the total weight 
of the structure in such a way that the flutter velocity and the 
divergence velocity are greater than or equal to some specified values. 

The problem can be stated mathematically as 

Minimize: w = f(P) 

subject to the restriction V(P) - ^0, U(P) I 

where P is the vector of variable design parameters P^ , V is the flutter 
velocity and U is the divergence velocity. and are the specified 
flutter velocity and divergence velocity constraints respectively. 

Since it is seldom possible or to optimize a structure such that 
V = Vq and U = simultaneously, the problem will be restated as follows: 

Minimize: w = f(P) (1) 

subject to the restriction <*> = (2) 

where <f> is either the flutter or the divergence velocity. 

An expression for the total weight of the structure can be formulated 
such that w is a linear function of the design parameters. 

Wo 

4 ~ ' 

where is the non-variable structural weight. Using the method of 
Lagrangian multipliers, the parameters Pj, P 2 , ... P^^ must be found 
which will cause the function 


( 3 ) 



to be an extremal subject to the condition 




Since F(P) must be an extremal 
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where w^. are a nstand derivatives of w with respect to P^. 
Let 'i’(P) b? approximated by the e..y.-ession 

(fie)- at2i.Pi 
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Substituting equation (6) into equation (5) yields 
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The parameters b^. and c^. are computed from the following relations 


- 

i)Pi, 




k- 

Ij + 


( 8 ) 

(9) 


"'+’1 ' d (^2 

where P. and T5- are values of P4 and respectively at point 1 in 
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the design space while P^ and a*"e values of P- and respectively 


at point 2 in the design space. 


From equations (7) and (8) it is seen that 
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and 

^ J = I 4-1 

where the values of a, b^ and c- will change with each redesign step. 
Substituting equations (6) and (7) into equation (4) yields the equation 

or '* 
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solving for a 



(14) 


Substituting equation (14) into equation (7) yields the following 
recu>"rence equation for P^: 

Pj = (m^- 

where b^ , and a are given by equations (10), (11) and (12). 




II. Qpliiii1i: j t1on Pro c edure 


Thu optimization procedure utilizes the recurrence equation (15) 
in a repative manner to compute new values nf P.. The procedure begins 
at soiiio arbitrarily cliosen point in the desigti space where and 
-f|V^ are computed. One gradient step is then made to find the values of 
f’ ■ and 0(|uations (10), (11), (12) and (15) arc used to compute 

* 'i.. 

new values of P.. Next I'.; , yj and are set equal to values of P. , 

1 ij ..r- ^ ’ 

j ■ 

; and , respectively and the design parameters P. are set equal to 

the values of computed from equation (15). New values of and 

are computed from values of P. and the process is repeated until 

all the values of approaches a small number e. When any of 

the design parameters become less than a minimum gage size, then Lpat 

parameter P. is set equal to the gage size with P. = P. - t.. When 
J ^ I ^ 2 

any value of P. is such that jP. - P. j < c then P. is set equal to 

1 1 2 Tv 

P- with P. = P. - c. Cor,.jtation is continued until |P,--P- I 

I ' 1 I ;> ^ 2 

approach some small number c. 

III. Application 

The method was applied to the problem of Ref. (1) which was a lifting 
surface supported by a simple box beam. The method appeared to be superior 
to tiiat of Ref. (1); however, the comparison was not conclusive since an 
error in the program for finding the stiffness matrix was discovered. 

The error will be corrected. The program used for finding the stiffness 
matrix was not the same as that used for Ref. (1). 

At this writing work is underway to apply the program to a fiber 
reinforced wing. The wing is straight with a unifonn cross-section and 
zero sweep angle. The fiber directions are fixed in four different 



directions 0”, 45°, -45° and 90° with respect to the span. There are 
eight lamina in the laminate which are arranged synmetrical ly about the 
midplane of the laminate. Each lamina has one fiber direction. The 
thickness of the lamina are the design parameters for the structure. 

A program for the analysis of the structure has been completed. 

IV. Conclusions 

The optimization method which was presented here is a workable 
method; however, its computational efficiency and usefulness have not 
been established by testing it on a structure with many design parameters. 
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The optimization procedure utilizes the recurrence equation (15) 

In a repatlve manner to compute n©rt values of P^. The procedure begins 
at son® arbitrarily chosen point in the design space P^^ where and 
are computed. One gradient step Is then made to find the values of 
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P^^ and ^2 equations (10), (11), (12) and (15) are used to compute 

new values of P.. Next P. , and ||i are set equal to values of P. , 
' '1 '2 

^2 snd , respectively and the design paraireters P^^ are set equal to 

the values of P. computed from equation (15). New values of ^2 and 
11^ are computed fran values of P| and the process is repeated until 
all the values of approaches a small number c. When any of 

the design parameters become less than a minimimi gage size, then that 

parameter P. is set equal to the gage size with P, » P. -e. When 
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P^ with P^ = P^ - e. Co:,.KUtation is continued until I 
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approach sane small nuntier e. 


III. Application 

The method was applied to the problem of Ref. (1) which was a lifting 
surface supported by a simple box beam. The method appeared to be superior 
to that of Ref. (1); hofifever, the conparison was not conclusive since an 
error in the program for finding the stiffness matrix was discovered. 

The error will be corrected. The program used for filling the stiffness 
matrix was not the same as that used for Ref. (1). 

At this writing work is underway to apply the program to a fiber 
r^ipforced wing. The wing is straight with a unifdra cr<»s*sectim and 
zero sweep angle. Tlie fiber directions are fixed in four different 
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directions 0“, 45®, -45® and 90® with respect to the span, There are 
eight lamina In the laminate which are arranged synmetrlcally about the 
midplane of the laminate. Each lamina has one fiber direction. The 
thickness of the lamina are the design paran^ters for the structure. 

A program for the analysis of the structure has been con^leted. 


IV. Conclusions 

The optimization method which was presented here Is a workable 
method; hmt/ever, Its computational efficiency and usefulness have not 
been established by testing It on a structure with nmny design parameters. 
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